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By interpreting the well-known, qualitative criteria for the existence of quantum error correction 
(QEC) codes by Knill and Lafiamme from a quantitative perspective, we propose a figure of merit 
for assessing a QEC scheme based on the average fidelity between codewords. This enables us to 
quantify the performance of a protocol as a whole, including errors beyond the correctable set. 
Various examples are calculated for photonic qubit codes dealing with the experimentally relevant 
case of photon loss, demonstrating the advantages of the new measure. In the context of continuous- 
variable QEC, our codeword-overlap measure can be used to reproduce, in a different way with 
no need for calculating entanglement measures, a previous result concerning the impossibility of 
improving transmission of Gaussian states through Gaussian channels via Gaussian operations alone. 
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I. INTRODUCTION 

Operating a device or scheme in the microscopic do- 
main places formidable demands on the purity and sta- 
bility of the diverse systems involved, augmented by 
fundamental limitations of Quantum Mechanics !]■ 
this context, ingenious mechanisms to avoid the unde- 
sirable effects of decoherence play a fundamental role in 
many areas of Quantum Information Processing (QIP), 
most notably, enabling large-scale, fault-tolerant quan- 
tum computation [2] and the communication of quantum 
bits across large distances [3] . Such mechanisms typically 
rely on two different solutions, namely, those provided by 
Quantum Error Correction (QEC) [1] and teleportation 
[5] combined with Entanglement Distillation or Entan- 
glement Purification Protocols (EPPs) 15]. 

As these techniques mature over a multitude of imple- 
mentations, with many advanced protocols already find- 
ing applications in real-world conditions, the necessity to 
define suitable measures to compare and rank different 
schemes becomes of uttermost importance. In particu- 
lar, the need arises to find quantities which can be easily 
computed for the widening range of experimentally acces- 
sible states currently used in QIP. Two main measures, 
in different variants, currently share widespread accep- 
tance: the entanglement degree is traditionally seen as 
the golden measure when defining the 'quantumness' of 
a channel, and is the benchmarking criterion of choice 
when evaluating EPPs and quantum memory devices [7] . 
Conversely, in the realm of QEC, fidelity measures, in 
particular, the worst-case fidelity or entanglement fi- 
delity [S], have historically been employed. The connec- 
tion between these approaches has been investigated in 
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different regimes [TO] , but it is not always a straight- 
forward one: Naively, one could expect codes which im- 
prove such fidelities to lead to higher safeguarded entan- 
glement, which is indeed true in certain cases; neverthe- 
less, there are examples [TTI - fTO] where a code will not 
safeguard any entanglement, up to and including caus- 
ing entanglement sudden death (ESD) [T3], while still 
improving input-output or entanglement fidelities under 
certain conditions. 

In this contribution, we propose a new measure based 
on the average fidelity between codewords of a given al- 
phabet, emphasizing that it satisfies certain desirable fea- 
tures. We present initially a motivating argument start- 
ing from well-known criteria for the existence of QEC 
codes, and proceed with examples comparing the entan- 
glement (as measured, in the case of qubits, by the con- 
currence |15| ) to our average codeword fidelity for a vari- 
ety of photonic schemes, with particular attention given 
to codes designed to protect a logical qubit against am- 
plitude damping [16] . In the spirit of [17] , we go beyond 
the analysis traditionally restricted to the correctable er- 
ror set and establish a comparison method for general 
schemes. However, here we impose no limitations to a 
code's distance or a priori assumptions on the operating 
regime [T^ (correction or detection modes ^9j). Finally, 
in the context of continuous-variable QEC with the logi- 
cal states being infinite-dimensional, we employ our new 
measure to re-obtain a known result I20j concerning the 
impossibility of performing QEC for Gaussian signals and 
channels when restricted to a Gaussian toolbox. 



II. ERROR CORRECTION CONDITIONS 

Consider a completely positive and trace-preserving 
map £ corresponding to a - possibly noisy - transmis- 
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sion/evolution channel. Its action on an arbitrary state 
Q is given by 

£{s) = Y^A^eAl, (1) 

fe 

dubbed the Kraus operator-sum representation pi] . 
Based on this channel decomposition, the Knill- 
Laflamme (K-L) criteria [52] establish that an error 
correction map TZ satisfying TZo £ = Id exists provided 
it holds that 

{Xi\AlAi\xj)=S{i-j)Xk,i, (2) 

where \xi) are the codewords from a given input alpha- 
bet, and the Xk^i define how different errors skew the 
codespace; in a simplified interpretation, it requires that 
different codewords remain orthogonal after the action of 
the channel, and that different errors must effect the same 
deformation across the input alphabet. Exact satisfiabil- 
ity of the above can only be achieved in certain scenarios 
[23]; at the same time, allowing for small deviations in the 
orthogonality and deformability requirements - or, alter- 
natively, ioT TZo£ to be close, but not necessarily equal, 
to the identity [24[ [25] - enables one to obtain more effi- 
cient codes [121 [25] . The discrepancies may be employed 
to bound the entanglement fidelity obtained with the use 
of the code [27]. 

Before proceeding, a few notes are in order. First, in 
achieving the codewords, an encoding step is implicitly 
assumed, and, for notation purposes, it is condensed to- 
gether with the transmission map into a single channel 
in the case of no encoding/direct transmission, one can 
simply regard the encoding as the identity map, Td. Fur- 
thermore, we note 6(i — j) corresponds to the more gen- 
eral case of a continuous alphabet; in the case of qubits 
(or other discrete- variable codewords), this should corre- 
spond to a Kroenecker delta, Sij. 

Motivated by the idea of employing a primarily quali- 
tative criteria as a potentially quantitative measure [28] , 
we will explore taking into account violations from the 
exact satisfiability of Eq. ^ (instead of neglecting such 
violations up to a certain order in the channel parameters 
[12j). Deviations from Eq. ^ can be broken down in two 
qualitatively different types: first, violations of 6{i — j), 
which lead different codewords to overlap, and thus re- 
duce the distinguishability between the input alphabet; 
second, departures from Xk^i, which in the non- violated 



^ Explicitly, this means that the initial signal Hilbert space is ex- 
tended, typically by adding a sufficient set of ancilla states and 
applying a global encoding unitary on the signal-ancillae system. 
The total system is then subject to a global channel transmission, 
which is usually composed of the original channel acting individ- 
ually and independently on the signal and ancilla subsystems. 
In the optical setting, enlarging the Hilbert space may mean ei- 
ther adding extra modes or allowing for higher photon numbers 
with the same number of modes. In the latter case, there are no 
additional channels that must act on auxiliary modes. 



case is strictly independent of the codewords, but in the 
violated case may affect different codewords in unequal 
manner (and hence may deform a superposition of code- 
words). Obviously this latter deformation can only take 
place when the prior effect is also present, but the con- 
verse is not necessary]^ We'll dub codes (i.e., encoding 
together with transmission channels) in which the en- 
tire alphabet is affected uniformly as "non-deformable" , 
whereas those with codeword-dependent skewness will be 
called "deformable" . 

Non-deformable codes will not preclude the distin- 
guishability to decrease, neither prevent different pairs 
from suffering varying deformations: the label only guar- 
antees, for any two orthogonal states, the overlap be- 
tween the original and the resulting states to be the same. 

In the case of experimentally relevant amplitude damp- 
ing channels, direct transmission of the |0) and |I) states 
(taken here as the occupation numbers of a bosonic mode 
like in so-called single-rail encoding) yields a deformable 
code: While the vacuum state is unaffected, the single- 
excitation 1 1) is taken to (1 - 7)|I)(I| -h 7|0)(0|. That 
the conjugate basis would be uniformly affected does not 
alter the classification of the code; in fact, for every de- 
formable code there exists at least one superposition in 
which the deformations are equal. However, considering 
the dual-rail encoding, |0)l |0)|1) and |1)l = |f)|0) - 
and assuming, of course, equal dampening in both rails 
(modes) - one finds again a non-deformable encoding, 
as any superposition of the codewords will result in an 
equal, global reduction of the length of the logical state 
vector. 

With the above considerations in mind, one could con- 
ceive employing the amount by which the criteria in 
Eq. ([2]) have been violated as a measure to rank different 
channels. We aim at a measure which, instead of classify- 
ing codes by the number and kind of errors it is capable 
of handling, should deliver information precisely about 
the operators the scheme is not capable of correcting, 
quantifying just how much it causes the input alphabet 
to skew and overlap. To this purpose, we construct the 
codeword overlap of a map by considering the output fi- 



^ More precisely, when taking each channel Kraus operator sep- 
arately like in Eq. of course each Kraus effect can af- 
fect different codewords differently, even though their overlap 
remains unaffected. An example is the amplitude damping 
channel acting on a single-rail qubit basis, for which we have 
(0|Ajylo|0> = Ao,o = 1, (1|aJAo|1> = Ao,o = 1-7, but also 
{0|Ajylo|l> = (see later, Eq. |4|). In other words, the length 
of vector |1} is reduced for any non-zero 7 and that of vector 
|0) remains unity, while they would perfectly maintain their or- 
thogonality. However, note that by taking into account a second 
Kraus effect, which together with the first Kraus effect forms 
a trace-preserving qubit channel, the deformation is indeed ac- 
companied by a reduction of the codeword distinguishability: 
(0|AqAi|1) = .^7. Our codeword overlap measure, as defined 
in Eq. ([sj for the whole trace-preserving channel, therefore cap- 
tures both deformation and non-orthogonality at the same time 
(see below). 
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delities obtained from a pair of orthogonal input states, 
and then averaging this quantity over the entire input 
alphabet: 



7CW 



dQ Ti 



/PQ 



(3) 



Here, pg and pg correspond to the outputs originat- 
ing from a pair of orthogonal input states \Q) and 
with \Q) _L \Q). The above loosely corresponds to an 
integration of Eq. ([2|, with two important distinctions: 
first, since the performance of a code should be basis- 
independent, we cover the surface of the Bloch sphere, 
taking each pair of diametrically opposed states as possi- 
ble codewords; second, instead of acting with the Kraus 
operators individually, we consider the full channel's ef- 
fect on the codewords, which enables it to be employed 
even in channels were a decomposition in form of Eq. ([T]) 
is not known. 

When considering encodings for spaces larger than 
qubits {e.g., the qudit codes proposed in [16]), one should 
per - doneform the integration considering opposing 
states in the surface of the corresponding Bloch hyper- 
sphere This, however, precludes a straightforward 
generalization for continuous-variable (CV) encodings; 
nevertheless, one may consider only a truncated set, i.e., 
a qudit encoding where d — )■ 00. 



III. PHOTON-LOSS QUBIT CODES 

Here we consider the ubiquitous binary logical basis, 
that is, with the information encoded in orthogonal states 
\0)l and The transmission of such states via opti- 

cal fibers corresponds generally to a lossy process, with 
the fiber absorbing or "losing" photons along its length. 
This corresponds to a Gaussian channel of particular rel- 
evance, dubbed an amplitude damping channel. It is 
characterized by a loss parameter, 7, and described, in 
the Kraus representation of Eq. ([T]), by an infinite sum 
{k = O...00) of operators of the form 

^'^ = E y^\/(l-7)"-V1n ~k){n\. (4) 

Diff'erent encodings have been designed to protect 
qubits against the errors caused by these operators, and 



^ The form above is of course not exactly convenient when consid- 
ering such encodings of higher dimensions. Already in the case of 
a qutrit (d = 3) this is easily seen: an arbitrary logical state has 
the general form a|0)ij -I- -I- c|2)i, and a possible codeword 

is lO)^, which has infinitely many states orthogonal to it (for 
instance, any state of the form -I- c\2) l with + = 1). 

Thus, in order to calculate the measure, for every |Q} state, one 
has to consider the average fidelity to each of the possible \Q) 
orthogonal states. 



are presented below to illustrate the usage of the new 
measure. The first code considered is the aforementioned 
dual-rail encoding, which provides detection capabilities 
only: 



|0)i = |01), |1)l = |10) 



(5) 



Whenever both modes are found to be in the vacuum 
state, the information has to be transmitted anew (error 
detection mode) or, after the decoding step, the output is 
replaced by the mixed qubit p = 5|0)(0| + (error 
correction mode, which will be considered here). 
The ubiquitous three-qubit repetition code, 



\0)l = 1000) , |1)l = |111) , 



(6) 



is capable of correcting any single bit-fiip error, and can 
also be employed against amplitude damping errors, as- 
suming photon loss to be a particular case of bit-flip, with 
a highly asymmetrical behaviour in which only one of the 
logical states is affected. 

Codes can also be constructed by exploring higher oc- 
cupations of the bosonic modes. For instance, in [16, , the 
authors develop the following encoding, capable of cor- 
recting the loss of up to one quanta to the environment: 



1 40) + 1 04) 
V2 



= |22) 



(7) 



Allowing for codes which do not exactly satisfy Eq. ^ 
enables one to achieve more economical encodings. One 
such was proposed by Leung et al. |12| . 



|0)l = 



10000) 



mil) 



V2 



\1)l = 



10011) + |1100) 

71 



, (8) 



which, assuming the damping constant 7 be kept low, is 
capable of correcting the loss of one excitation to 0(7^). 

Finally, the usage of infinite-dimensional carriers is one 
that offers numerous advantages in the realm of QIP (SUl 
ISlj . One may thus consider the encoding of qubits on 
coherent states, which are innately non-orthogonal: 



a 



ID 



(9) 



Here, the encoding can be seen as a channel which ir- 
reversibly reduces the distinguishability; in this case, a 
state is no longer orthogonal to the state found on the di- 
ametrically opposing point of the Bloch sphere, with the 
exception of those lying in the equator. Error correction 
codes have also been developed in this regime, such as 
the protocol in i32 , which effects, up to a normalization. 
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a) 
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(10) 



capable of correcting ['^^^J amplitude-damping induced 
phase-flip errors. We note, however, that not only the 
transmissivity of the channel contributes to this scheme's 
performance: also the size of the superposition, which 
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FIG. 1. (Color online) Codeword overlap (top) and concur- 
rence (bottom), as a function of the damping parameter 7, 
for different codes: dual-rail, Eq. (5) (blue), three-qubit rep- 
etition, Eq. (6) (red), bosonic, Eq. (7) (green, dashed), and 
four-qubit, Eq. (8) (black, dotted) [^ . 



regulates the amount of overlap in the input alphabet, 
affects the transmission characteristics (see e.g. [3.3^ 3|]). 

We proceed by calculating the average codeword over- 
lap for each of the above encodings and comparing it to 
a different figure of merit, namely, the amount of entan- 
glement preserved after employing the protocol at hand 
to transmit one half of a two-qubit maximally entangled 
state (further details are provided in Appendix A). The 
measures are computed as a function of the channel loss 
parameter 7, in the case of the discrete-variable encod- 
ings, or for a fixed channel transmissivity, as a function 
of the coherent-state superposition size \a\. 

The results are presented in Figs. 1 and 2. In the for- 
mer, we observe that the overlaps are zero when the chan- 
nel approaches perfect transmission: as expected, the 
orthogonal input codewords remain orthogonal; equally, 
the entanglement remains maximum. As the losses in- 
crease, the overlaps grow; at the same time, the entangle- 
ment diminishes. Alone, this fact is highly unsurprising; 
however, we remark that also the ordering established 
through one measure is reflected on the other. That is to 
say, the code which offers the best performance, in terms 
of codeword fidelities, in a given regime, also returns 
the highest safeguarded entanglement; as the parame- 
ters vary, the relative ordering changes as well, and this 
is observed in both figures of merit. However, the exact 
point in which a change of ordering occurs is not strictly 
always the same: between certain encodings, the chan- 
nel parameters in which a crossing occurs may differ by 
small amounts. This is particularly noticeable when ex- 
amining the three-qubit repetition codes, which is highly 
deformable and suggests that variations may be related 



FIG. 2. (Color online) Codeword overlap (top) and con- 
currence (bottom) for direct transmission Eq. (9) and the 
coherent-state code Eq. (10), after transmission through a 
7 = 0.32 channel, as a function of the coherent-state super- 
position size \a\. Encodings with 3 (red), 5 (green, dashed), 
11 (black, dotted), and 51 (grey, dashed) qubits are depicted 
in addition to the direct transmission case (blue) [SJ . 



to the codeword-dependent skewness. 

All of the above suggests a strong link between the 
overlapping properties of the output alphabet, and the 
amount of entanglement capable of being safeguarded by 
means of a given scheme. Nevertheless, we note that no 
a priori reason exists to suggest that the crossing points 
should be precisely the same: after all, we are dealing 
with measures of apparently different character. 

A further observation is in order: while certain encod- 
ings lead to entanglement sudden death [TTJ [131 HI! > 
equivalent catastrophic breakdown in terms of overlaps 
(average codeword fidelity equalling unity) does not hap- 
pen in any encoding until the channel becomes fully lossy. 
This would suggest that the new measure is capable of 
portraying certain characteristics which would be other- 
wise lost in an analysis solely based on the entanglement. 

In Fig. 2, we observe the behavior of different coherent- 
state encodings as a function of the superposition size 
|a|. Again, a good agreement is found between the or- 
dering obtained through concurrence and codeword over- 
laps. Notable is the difference in the case of direct trans- 
mission: while it is possible to obtain a maximally entan- 
gled state in the limit of a — >■ 0, the alphabet as a whole 
becomes indistinguishable and thus impractical for en- 
coding purposes - a feature reflected by the codeword 
overlap measure. 
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IV. GAUSSIAN ERROR CORRECTION NO-GO 

Optical modes of the electromagnetic field, whose 
quadratures satisfy the canonical commutation relations, 
provide a natural testbed for a wide range of quantum 
information concepts |35j . Here, logical information is 
encoded by means of a truly continuous logical alphabet; 
for instance, the codewords defined by the position eigen- 
states {la;)^}. In this scenario, the so-called Gaussian 
operations - defined as those that map Gaussian states 
into Gaussian states - are of great relevance due to the 
ease in which they can be implemented with current ex- 
perimental resources. 

Nevertheless, the capabilities of the Gaussian set are 
not without restrictions. Given the fragility of the quan- 
tum resources in face of ubiquitous decoherence mech- 
anisms, a significant limitation is that such operations 
are incapable of distilling higher entanglement from less 
entangled Gaussian states |5SH55] . or, in close relation, 
unable to protect Gaussian states from the widespread 
class of Gaussian errors [20] . 

Interestingly, Gaussian transformations alone sufRce to 
suppress and correct non-Gaussian errors acting on arbi- 
trary states, but in particular, the Gaussian transforma- 
tions also suffice when the non-Gaussian errors act upon 
Gaussian input states [55HI^ ; equally worth noting is the 
fact that a Gaussian error channel, when acting upon spe- 
cific non-Gaussian input states, can exhibit stochastic, 
non-Gaussian behaviour |43j . however, in this case Gaus- 
sian encoding and decoding procedures alone appear to 
be incapable of correcting such errors, and non-Gaussian 
operations should be accounted for [32] . 

In this section, we'll employ the codeword overlap mea- 
sure in order to explore what kind of statements it allows 
us to make in the context of continuous- variable QEC, in 
particular, in relation to the aforementioned No-Go the- 
orem for all-Gaussian QEC [5D]. Most importantly, such 
an approach enables us to obtain fundamental insights 
in a rather distinct way, without the need for calculating 
an entanglement measure (such as the logarithmic nega- 
tivity employed in Ref. [2D])- We note that the quantity 
as developed in Eq. (|3| is based on averaging the fidelity 
between the channel-output states which originate from 
two orthogonal input states. However, here, in order to 
include the possibility of non-orthogonal input states as 
it is typically the case for a Gaussian alphabet, we'll con- 
sider two general states: if the output overlap between 
two arbitrary states cannot be reduced through QEC, 
then, in particular, it also holds that the overlap between 
two orthogonal states (after their channel transmission) 
will not decrease, and integrating over all possible pairs 
of states will furthermore not result in a lower quantity, 
thus establishing the desired result. 

The fidelity between two states pi and p2 is given, in 
the characteristic function formalism, by 



F 



1 



For Gaussian states with zero mean, the characteristic 

1 T 

function is given by Xpi = e~^^ and we can re- write 
the above expression as 



F : 



1 



which evaluates to 



(12) 



(13) 



with A — det (tri + (T2) and 6 — (detai — 1) (det (72 — 1). 
Now, the action of a Gaussian channel on the level of the 
covariance matrices is given by 7 — >■ MjM"^ + N (see 
Appendix B for details). Taking this transformation in 
acount, the resulting fidelity is 



F' = 



2 



(14) 



(11) 



where A' = det {MaiM'^ + Ma2M'^ + 2N) , and 
(5' = (det(McriM^ + N)-l) {det{Ma2M'^ + iV) - l). 

Now, to establish a relationship between F and F' , 
we must evaluate the resulting effect from M and N in 
the above (for detailed arguments see Appendix C). We 
are faced with three relevant cases: (i) |detM| = 1, (ii) 
I detM| > 1 and (ui) | detM| > 1. The first case is triv- 
ially evaluated: when det = 0, this corresponds to a 
Gaussian unitary, which, as expected, yields exactly the 
original fidelity, i.e., F' = F. For detA^ > 0, represent- 
ing the addition of classical (thermal) noise, one observes 
the difference between the square roots in Eq. ( 14 ) to di- 
minish, and thus the fidelity to increase. The second case 
also finds a straightforward solution based on the same 
argument. The third case, however, contemplates chan- 
nels which do not necessarily induce a spreading of the 
Gaussian state, and requires a more careful analysis. In 
this case, the action of the map results in a "contraction" 
towards a common state. This is clearly exemplified by 
the prototypical amplitude damping channel, using two 
displaced thermal states as inputs: as the loss parameter 
7 increases, the states are gradually attracted towards 
the vacuum state; however, the increased purity (and 
correspondingly reduced values of det ai) plays no role in 
diminishing the overlaps. 

We note furthermore that displacements {i.e., shifts 
in the first-order moments) bear no effect in the above, 
simply rescaling F and F' by a fixed amount. 

One can then conclude that, except when £ corre- 
sponds to a symplectic operation, either the additional 
noise will cause a spreading of the Gaussian state, or the 
states will be contracted towards a common state. In any 
case, the fidelity to any other state subject to the same 
action is bound to increase. When the decoding opera- 
tion is equally symplectic, one can at best re-obtain the 
original state (through TZ — S^^), assuming, of course, 
that the encoding operation already had the best choice 
of unitary operations. 
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In the above discussion, eventually we considered in- 
put states of a Gaussian nature, thus reproducing the 
known No-Go result for all-Gaussian QEC [5D] in terms of 
our fidelity-based codeword criterion and independently 
of entanglement measures. The treatment of Gaussian 
codeword pairs above is analogous to our earlier treat- 
ment for the qubit codes, and the argumentation follows 
through depending on the assumption of a Gaussian in- 
put alphabet (see Appendix C). 

V. CONCLUSIONS AND OUTLOOK 

We have defined a quantitative measure based indi- 
rectly on the amount by which the well-known quantum 
error correction criteria by Knill and Laflamme are vio- 
lated, proposing the use of averaged codeword fidelities as 
a figure of merit to evaluate different quantum protocols. 
By quantifying how distinguishable originally orthogonal 
inputs emerge from an error channel, the measure finds a 
natural interpretation as a translation, into the quantum 
regime, of the classical coding theory notion of confus- 
ability of an alphabet [15] . 

While initially harder to compute, by requiring the 
evaluation of the fidelity between two generally mixed 
states, the codeword overlap measure does not require 
the optimization (minimization) necessary to obtain the 
worst-case fidelity. At the same time, it accurately de- 
picts the alphabet's distinguishability behaviour which 
fails to be portrayed by means of the entanglement fi- 
delity or conventional entanglement measures. We em- 
phasize that properly describing this distinguishability, 
both in an ideal and in a noisy or lossy quantum channel, 
becomes particularly important when the ideal (quan- 
tum) codewords are already non-orthogonal, as it is of- 
ten the case when the quantum information carriers are 
continuous- variable oscillator states. 

The employment of the codeword overlap reflects qual- 
itatively the behaviour found when quantifying the per- 
formance by means of other figures of merit, in particular 
here the concurrence. In other words, whenever the vari- 
ation of a certain channel parameter causes one scheme to 



improve (reduce) the average codeword fidelity in com- 
parison to another scheme, one also observes that the 
first scheme will result in an improved figure in the safe- 
guarded entanglement. This matching was also found to 
be quantitatively exact for certain choices of codes, how- 
ever, the precise reasons for such behaviour are hitherto 
unknown. 

We have also employed the new measure to reaffirm 
the impossibility of improving the transmission of Gaus- 
sian states, subject to Gaussian noisy channels, through 
Gaussian operations alone. In this case, the No-Go re- 
sult was found independently of an entanglement mea- 
sure, solely based upon the overlap of Gaussian codeword 
pairs. 

Finally, since the measure is, in principle, accessible 
and computable in those instances where entanglement 
is found hard to quantify (i.e., infinite-dimensional non- 
Gaussian states), we expect these results to shed light in 
other aspects of (photonic) error correction in particular 
and quantum information in general. 
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Appendix A: Calculating the concurrence and 
average codeword fidelity 



In obtaining the codeword overlap curves depicted in 
Figs. (1) and (2), the following procedure is used: 

For each encoding, an input qubit \Q) is prepared in 
the state cos ^|0)i + e*^ sin concomitantly, a sec- 



ond state is prepared, \Q) = sin^|0)_L 



cos ■ 



ensuring that {Q\Q) — 0. Each of those states is then 
subject to the transmission channel characterized by the 
operators in Eq. Followed by each code's respective 
decoding procedure, the following outputs are found: 
Direct transmission results in 



PQ , direct — 

1 / 1 + 7 + cos W 
and 



7COSW e ^/l — ^svnw 
7sinii; (7 — l)(cos w — 1) 



(Al) 



P Q, direct 



- cosw + (1 — 7) cos w — e 
-e'^v^l — 7sinu' —(7 



-ie 



v^l — 7 sin w 
l)(cosu' + 1) 

(A2) 



To the dual-rail encoding corresponds 

PQ. dual —rail — 

1 (1 + (-1 +7)cosw) e"*^(-l -|-7)sinu; 
2\ e'^(— 1 + 7) sin w (1 + cos w — 7 cos w) 

and 

Pq, dual — rail 

1 / (1 + cosw — 7COSW7) — e^*^(— 1 + 7) sin w 

2 V -e'''(-l + 7)sinw; (1 + (-1 + 7) cos w) 



(A3) 



(A4) 



The 3-qubit code produces 



and 



_ 1 / 1 + (3- 2p)p2 + (p- 1)2(1 + 2p)cosio e-*^(l -p)3/2sinw ; 
PQ,3-quMt- e^«(l-p)3/2sin«; (p - 1)^(2 + 4p) sin^ f ' 



^ 1 /l + (3-2p)p2_(p -1)2(1 + 2p) cos u; -e-^»(l - p)^/^ sinu; . 
PQ,3-quMt 2\ -e'''(l -p)3/2sinw (p - 1)2(1 + 2p)(l + cosw) / ■ ^ ' 
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For the bosonic encoding, one finds 



PQ^bosonic 



and 



1 1 - (7- 1)^(1 + 37) cos u; -e-*^(7- 1)3(1 + 37)sinw7 
~e^^(7 - 1)3(1 + 37) sin u; 1 + (7 - 1)^(1 + 37) cos w; 



_ 1 / 1 + (7- 1)3(1 + 37) cos w e-*^(7- 1)3(1 + 37) Sim 



PQ,bosonic 2 e'^(7 - 1)3(1 + 37) sinw 1 - (7 - 1)3(1 + 37) cos w 
Finally, for the "approximate" encoding, the outputs are 

PQ , approximate — 

i (1 + 7^(27 - 1) + (7 - 1)2(1 + 27) cos w) ie*^ (7^ - 73 + er"^'^ (2 + 7^(37 - 5))) sin w 
y-^e (-^2 _ ^3 _^ g2*e (-2 + ^2(3^ _ 5)-)) gj^^^ 1 (l + ^2 _ 273 - (7 - 1)2 (1 + 27) cosw) 



(A7) 



(A8) 



(A9) 



and 



^approximate 



(l + 7'( 
((-1 



4 + 27) - (-1 + 7)2(1 + 27) cosiw) 
' 7)72 + e2'« (-2 + (5 - 37)72)) sin w 



s-'«(-l + 
1(1 + 7^ 



7) (2 + 7 I 
- 273 + (^ 



2+ (-3 + e2''») 7)) sinw 
1 + 7)2(1 + 27) cosu;) 



(AlO) 



r 



With the above, Eq. ([3]) is then computed by means 
of a numeric integration procedure over w and 9. At 
each value of 7, approximately 1000 states are employed 
to produce the average, although for less deformable en- 
codings a smaller sample already resulted in an adequate 
agreement with the asymptotic behaviour. 

For evaluating the entanglement-safeguarding capa- 
bilities of a scheme, the maximally entangled state 
(|0)|0)l + |l)|l)i) is employed, transmitting the en- 
coded mode through the lossy channel, afterwards fol- 
lowed by the decoding procedures. We note that the 
decoding operations reduce the encoded state back to a 
qubit subspace, therefore allowing the entanglement to 
be computed on a 2 x 2 Hilbert space. 

For the entanglement analysis developed in the text, 
we consider Woofers' concurrence, which is obtained 
through 



/aI} . (All) 



C — max{0, v Ai 



Here, are the eigenvalues, in decreasing order, of pp, 
where p = {uy^i ® <Jy^2)p* {<yy,i ® o'j/,2), and Uy^i is the 
Pauli Y operator in the i-th mode. 

In the case of coherent-state encodings evaluated in 
Fig. (2), the input state for both direct transmission and 
the different encodings is 



\Q) 



1 



./ma) 



{y/w\ -a) + e'Vl - w\a)) 



(A12) 



with 0>u'>l,0>6'>7r, and the normalization con- 
stant iV (a) = 1 + 2 cos Oyjwil - w)e-2|"l' (a6* +a*b), re- 
quired due to the non-orthogonal nature of the coherent- 
state alphabet. When considering the error correction 
protocol [351 , the deterministic operation of the scheme is 



considered; in other words, the transformation in Eq. ( 10 ) 



is only obtained for large superposition sizes, with smaller 
\a\ incurring an erroneous component (see ref. |34j for de- 
tails) . This induces in a double-tradeoff between the su- 
perposition size and the figure of merit being considered: 
while on the one hand, larger "cat states" are desired - 
not only due to their inherenty larger distinguishability, 
but also for the operation of the non-Gaussian Hadamard 
gates found in the protocol - on the other hand, this ef- 
fects an increased channel-induced phase-flip probability. 

For computing the concurrence in the later case, two- 
mode maximally entangled states are employed. 



1$-) 



1 



V2 - 2e-4|"l' 



|a,a) - I - a,-a)) . (A13) 



Again, the first mode is kept while the second is sent 
through the error-correcting scheme. Finally, for the 
above calculations, it helps to express the density ma- 
trices in terms of an orthogonal basis jfa)}, such 
that 



\a) = Pc,\Ua) + l^c,\Va) 
- a) = Pc,\Ua) - yc,\Va) 



(AM) 



with 



Pa 



and Vn = \ — 
\ 2 J \ 2 

Appendix B: Gaussian Formalism 
1. Gaussian States 



Here we briefly review the Gaussian formalism, adopt- 
ing the convention from j36j . 
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Given the Weyl operators 
W{x) = 



(Bl) 



,Xn,Pn) , with the 



where x € IR^" and R ^ {Xi,Pi,. 
commutator relations [X^, P^] = i6jk, a Gaussian state 
p is defined as having a Gaussian characteristic function 
Xpi^) — tT[pW{x)]. An equivalent definition can be given 
in terms of the state's Wigner function, which can be 
obtained from the characteristic function by a Fourier 
transform, and is also Gaussian-shaped. One defines the 
first and second order moments d and 7, respectively the 
displacement vector and the covariance matrix (CM), by 
di = (xi) and jij = {xiXj +XjXi) — 2didj. In most cases, 
d can be set to zero without loss of generality. The CM 
satisfies 7 = 7^ > iJ„, where we define the symplectic 
matrix 



fc=i 



-1 

1 



(B2) 



The displacement vector and covariance matrix com- 
pletely determine the state p, whose density operator can 
be written as 

p = 7r-"/ da;e-3^^''^+*'*''^iy(x) . (B3) 

Of special relevance is the maximally entangled state 
(MES), corresponding to an infinitely squeezed two-mode 
squeezed state (TMSS), with CM 



lim 7(r) 

r— f 00 



(B4) 



where = coshrl and Cr = sinhrA are both 2n x 2n 
matrices, and 



A = diag(l, -1,1,-1,..., 1,-1), 



2. Gaussian Operations 



(B5) 



A Gaussian channel is defined as a map £ taking Gaus- 
sian states into Gaussian states, cf. p' — £{p). Following 
the Choi-Jamiolkowski isomorphism between completely 
positive maps and positive operators ^46] , to every Gaus- 
sian map £ there corresponds an operator E, 



E12 = lim (f®l) (10)12(01) , 



(B6) 



this equation allowing us to re-interpret the transmission 
of an arbitrary state p through f as a teleportation using 
E12 as the entangled resource state, i.e., 



£{p) CC tr2[Ef^P2] = tr23(^12P3l0)23('/'l). 



(B7) 



One should note that, since £ maps Gaussian states into 
Gaussian states, and |0) in Eq. (B6) can be taken as the 



limit of a Gaussian state, E must itself correspond to a 



Gaussian operator which, similarly to (B3), can also be 
written as 



E = 



dx ( 



"Fx+iD'^x-C 



W{x), 



(B8) 



with appropriately-defined CM F, displacement vector 
D a nd a normalization constant C. Now, by empl oyin g 
Eq. ( B7 1 and replacing E12 with the operator in Eq. ( |B8| , 
one can obtain the action of f on a general state, ie, 
obtain 7' and d' in. £ : p^^d — ^ P-y^d' from F and 

We now consider the set of operations which can be 
implemented by augmenting our system with additional 
(Gaussian) ancillary states, performing Gaussian unitary 
operations over the whole combined system and discard- 
ing (tracing over) the ancillas, thus obtaining the class 
of Gaussian completely positive trace-preserving (CPTP) 
maps. The action of such a map on a state with CM 7 
is given by 



7 MjM^ + N . 



(B9) 



with M real and A^ > real and symmetric. We must 
also have 



detA^> (detM-1)^ 



(BIO) 



lest the complete-positivity requirement be violated. The 
Gaussian operator corresponding to this operation has 
the CAI 



lim 

r— f C30 



M'^ArM + N M^Cr 



(Bll) 



The solutions to Eq. (BIO) with A = (adding no ex- 
tra noise) and det M — 1 (preserving the sum of areas) 
are dubbed symplectic transformations. In the quantum- 
optical context, these are exemplified by unitaries such as 
squeezers, phase shifters, or lossless beam-splitters. One 
can easily verify those to be the channels with minimal 
entanglement degradation [5D] , or equally, those preserv- 
ing the newly-developed codeword overlap measure. 

Finally, we note that the identity map, I, is obtained 
in Eq. by taking M = 1 and A' = (being thus a 



symplectic operation), and reduces the expression (Bll) 
to t he C M of the maximally entangled (r — )• 00) TMSS, 
Eq. (IB4|. 



Appendix C: Fidelity between Gaussian states 

The effect of a Gaussian channel on the level of the 
covariance matrices is given by 7 — > MjM'^ + N. In 
order to facilitate the treatment of such expression when 
evaluating fidelities, we note that the above Gaussian 



* The interested reader should follow sections II and III in ref. |36| 
for a complete discussion of this procedure. 
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channel can be parametrized by means of a singular value 
decomposition. This gives rise to an equivalent channel 
characterised by the matrices 



M' = SVMU and 
N' = SVNV'^S 



(CI) 



Since U, V and S are all symplectic, the overlaps between 
the original and transformed channel remain unaltered. 
Now, detM' = detM (and equally for N' and N); fur- 
thermore, without loss of generality we can choose M 
as proportional to the identity, i.e., M oc r/l. Doing so 
greatly simplifies the expression for the fidelity after the 
Gaussian operations, since one can then basically employ 
the original determinants, up to a scaling factor. 

Now, to establish a relationship between F and F' , 
we must evaluate the resulting effect from A I and N in 
the above. We are faced with three relevant cases: (i) 
|detM| = 1, (ii) |detM| > 1 and (iii) |detAf| > 1. 

The first case, |detM| — 1, with detA^ = corre- 
sponds to a Gaussian unitary and is trivially evaluated, 
yielding, as expected, exactly the original fidelity 

For detA^ > 0, representing the addition of clas- 
sical (thermal) noise, and equivalently for the second 
case, |detAf| > 1, with help of the parametrization in 
Eq. ( |C1[ ), one trivially observes the difference between 
the square roots in Eq. ( 14 1 to diminish, and thus the 



fidelity to increase. The later holds for quantum-limited 
maps 47J where N bounds Eq. ( B10[ ); and evidently for 
those cases where a second, (classical) noise channel fol- 
lows. 

The third and final case, however, involves a more elab- 
orate analysis. Using the expression for the symplectic 
transformation Eq. (CI) in Eq. (14), expressing the rel- 



evant quantities in terms of the parameter rj governing 
M and N, and employing th e phy sicality constraints for 
Gaussian channels (see Eq. (BIO) in Appendix B), one 



finds, after long but otherwise straightforward calcula- 
tions, that 



F- F' >0 



(C2) 



has an empty solution set; and thus, the fidelity is indeed 
bound to increase. 

Finally, in order to fully appreciate that the treat- 
ment with codeword pairs of Gaussian states developed 
in section IV is analogous to our earlier treatment for 
the qubit codes in section III, and also that the ar- 
gument follows through depending on the assumption 
of a Gaussian input alphabet, notice the following. A 
pair of continuous-variable codewords may be written as 



= J dxtpi{x)\x)L and \ip2) = J dx tp2ix) \x) l , sim- 
ilar to our \Q) and \Q) states for qubits, but with the 
general Bloch-sphere parameters (App. A) this time re- 
placed by two general wave functions V'i(^) 9'iid ip2ix); 
the states and \tfj2) are also no longer restricted 

to be orthogonal. Now assuming that the logical states 
I a;) i that span the continuous codespace are the result of 
an arbitrary, unitary Gaussian encoding operation act- 
ing upon a one- mode basis state \x) together with an 
arbitrary Gaussian multi-mode ancilla state | Gaussian), 
\x)l = [^Gaussian {\x) jGaussiau)), we obtain for the two 
input states. 



= / dx4,,{x)\x)L (C3) 

dxip^ix) J/caussian {\x) (g) |Gaussian)) 

= f/caussian (^J dx ^pi{x)\x) (K" | Gaussiau) 

with i — 1,2. Under the assumption of two Gaussian 
wave functions ipi{x) and ip2{x), it is guaranteed that 
the resulting state above is again a Gaussian state. How- 
ever, for two arbitrary wave functions tpi{x) and ip2ix), 
the resulting states can be non-Gaussian, and the aver- 
aging for obtaining our overlap measure would have to 
include non-Gaussian codewords as well (corresponding 
to the more general scenario of arbitrary signal states 
subject to Gaussian encoding/decoding operations and 
Gaussian error channels). In this latter case, the above 
analysis for the Gaussian overlaps would no longer suf- 
fice. Therefore, we must assume that ipi^x) and ip2{x) 
represent any pair of Gaussian wave functions, and the 
above arguments follow through. Note that this conclu- 
sion is robust against a change of the basis of the original 
signal state, |a;) — > \n) (with \n) being the photon number 
basis), in which case 

|^») =^Q,nki)L (G4) 

n 

= ^ Ci,« ^/Gaussian (|«) |Gaussian)) 



\n) (g) jGaussiau) ] 



and where f/Gaussian Hn) <Si [Gaussian)) is generally a non- 
Gaussian state, but C/caussian (Z]n Ci,„ \n) (g) |Gaussian)) 
is a Gaussian state, provided that j-f/'i) = X)n '^i,n k*) is a 
Gaussian state too. 



